The Higgs resonance in fermionic pairing 
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The Higgs boson in fermionic condensates with the BCS pairing interaction describes the dynamics 
of the pairing amplitude. I show that the existence and properties of this mode are sensitive to the 
energy dispersion of the interaction. Specifically, when the pairing is suppressed at the Fermi 
level, the Higgs mode may become unphysical (virtual) state or a resonance with finite lifetime, 
depending on the details of interaction. Conversely, the Higgs mode is discrete for the pairing 
interaction enhanced at the Fermi level. This work illustrates conceptual difficulties associated with 
introducing collective variables in the many-body pairing dynamics. 



The low-energy properties of many-body systems with 
spontaneously broken symmetries are encoded in the col- 
lective modes associated with the structure of the order 
parameter. The dynamics is usually dominated by the 
gapless Goldstone modes that restore the symmetry [lj, 
which makes them easily observable in variety of systems. 
In contrast, the massive Higgs modes resulting from 
the amplitude variations of the order parameter are hard 
to detect. Moreover, they do not necessarily exist when 
the order parameter is generated dynamically due to col- 
lective effects. Understanding how these modes emerge 
in a microscopic framework is essential for constructing 
field-theoretical descriptions operating with a partial set 
of dynamical variables. 

This peculiarity of the Higgs mode becomes evident in 
the BCS model of superconductivity [3j], where pairing of 
fcrmions with opposite spins and momenta is described 
by a complex-valued order parameter. Several collective 
modes distinct from the Higgs mode have been consid- 
ered in superconductors: the excitonic states 0,0], the 
high-temperature phase modes 0] , and also a phonon gap 
mode observed in materials with coexisting superconduc- 
tivity and charge-density waves @, At the same time, 
no direct signatures of the intrinsic amplitude mode have 
been found so far. The main difficulty has to do with the 
structure of quasi-particle excitations in the system. 

The role of quasi-particle states is clarified in the re- 
sponse to a small pairing perturbation. The induced dy- 
namics can be described as the dephasing of individual 
quasi-particle states excited in the vicinity of the Fermi 
energy, with the rate of dephasing proportional to the 
excitation energy. As a result one finds a power-law de- 
cay of the perturbation Equivalently, this behavior 
can be linked to the square-root singularity of the density 
of states at the gap energy Q . Recent extension of the 
analysis to non- linear regimes (nil . 11 1 has also demon- 
strated the absence of the intrinsic collective excitations 
of the pairing amplitude. One of the implications of this 
behavior is the lack of an effective description of the pair- 
ing dynamics such as time-dependent Ginzburg-Landau 
theory in gapped superconductors [12, EH. However, the 
results of the analysis rely on the key assumption of con- 
stant interaction between fermions and change qualita- 
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FIG. 1: The Higgs mode energy sketched as a function of 
dispersion: the mode is discrete within the spectral gap for 
positive dispersion (regime A); the mode is non-physical (vir- 
tual state) for a smooth negative dispersion (regime B); and it 
becomes a resonance inside the continuum of the BCS quasi- 
particles for sharp and negative dispersion (regime C). 

tively otherwise. 

In this work, I discuss three different regimes of the 
pairing dynamics that appear due to the energy disper- 
sion of fermionic interaction in the BCS model. The main 
results of this work are summarized in Fig.[TJ In particu- 
lar, the energy of the discrete amplitude mode lies inside 
the spectral gap (regime A) resulting in non-decaying dy- 
namics when the pairing is enhanced at the Fermi energy 
(positive dispersion) . Upon smooth reduction of the dis- 
persion, the energy of the mode gradually approaches the 
continuum of quasi-particle excitations eventually merg- 
ing with the edge of the spectrum at constant interac- 
tion. In this case, and also for a smooth suppression 
of the pairing (negative dispersion) the collective mode 
becomes non-physical (virtual state), which is reflected 
in a power-law decay of the dynamics (regime B). The 
mode becomes a resonance (damped oscillation) within 
the quasi-particle continuum with a finite lifetime only 
when the energy-scale of suppression is small compared to 
the equilibrium gap at the Fermi energy (regime C). The 
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lifetime of the Higgs mode is small compared to the quasi- 
particle relaxation time, making the dynamics analogous 
to the Landau damping 14] in collisionless plasma. 
Our system is described by the BCS Hamiltonian 



n = Yl e P a P<r a P^ _ Vq«pT a -p| a ql a -qT 



(1) 



pa 



pq 



where a p<7 is the creation operator of a spin- 1/2 fermion 
with momentum state p and spin state a =t,l, £p = 
p 2 /(2m) — /i is the energy of free fermions with respect 
to the Fermi energy /i. The summation is conventionally 
performed over the states within the band |e p | < ll>b of 
finite width related to e.g. the Debye energy. 

The energy dispersion of interaction A pq in Eq.JT]) 
leads to the energy-dependence of the pairing gap. The 
analytical treatment is simplified in the case when it 
is taken in a factorized form A pq = A/ p / q , such that 
/| p | =PF = 1. Diagonalizing the BCS Hamiltonian J!}, 
one obtains the quasi-particle excitation spectrum E p = 
(e p + Ap) 1 / 2 characterized by the pairing amplitude 



(2) 



where the energy dependence is provided by f p , and A F 
is the equilibrium value of the pairing gap at the Fermi 
energy. The self-consistency condition takes on the form 
2/Ao = J2 P fp/E P and defines Aj?asa function of f p and 
the coupling constant Ao . The factorized form of energy- 
dependent interaction simplifies the analysis. Nonethe- 
less, the qualitative results obtained in this work also 
apply to a generic form of interaction. 

We analyze the dynamics induced by a small pertur- 
bation of pairing interaction, 6\(t) — X(t) — Xq. It is 
characterized by the expectation value of the pairing op- 
erator 



A = A ^ / P a pT aI 



p 



pI' 



(3) 



calculated with respect to time-dependent many-body 
state. The problem is solved by employing the time- 
dependent mean-field approach [15| valid at time scales 
smaller than the quasi-particle relaxation time. The lin- 
ear response of the pairing amplitude <5A p = f p SA is 
given by 



SA{t) = 



dr([A(t),V(r)}), 



(4) 



where the operators evolve according to the unperturbed 
Hamiltonian |T|), and the average (...) is taken with re- 
spect to its ground state. The perturbation operator 



V(t) 



A F /8\{t) 
Ao V Ao 



6A{t)\ 
A F J 



A^+h.c. (5) 



includes a term cx 5A(t) obtained from the averaging of 
the interaction energy within the mean-field approxima- 
tion. 



The solution of Eq. ([4]) simplifies if one assumes the 
particle-hole symmetry. This condition implies that the 
energy-dependent pairing gap is symmetric with respect 
to the Fermi energy. In this case, the phase and the 
amplitude modes of A decouple and can be analyzed 
separately. Concentrating on the amplitude mode 5 A' 
excited by the perturbation ([5]), we obtain the solution 
of the inhomogeneous integral equation by the Fourier 
transformation: 



SA'(lu) = [G(to) - 1] ^5X(w), 
Ao 



(0) 



where G(lo) is a retarded spectral function describing the 
many-body response. By the relation G -1 (a->) = 1 + K(u) 
it is expressed through a correlation function K{ui) that 
reads in real time as 



iK(t) 



Ao 



([A(i),At(0)]) + ([A(t),A(0)])). (7) 



Calculating K(t) we derive the expression for the spectral 
function 



G~ 1 M=. 



de 



u 2 /4- 



A c 



A F J E e (EI-u> 2 /4Y 



(8) 



where g = i/\q is dimensionless coupling, and 
the particle-hole symmetry is enforced by energy- 
independent density of states v. The states in © are 
labeled by energy e instead of momentum p, and accord- 
ing to the usual rule u> — > u + iO for the retarded func- 
tion. It is worth noting that Eq. ([8j) is in agreement with 
Rcf. 0. 

The analytic structure of G{oj) in the complex w-plane 
defines the time-dependence of the pairing amplitude 
8A{t). In particular, we are interested in locating the 
poles oj a = ujq — iT of G{uj) cx (u> — w a )~ 1 , or equiva- 
lently the simple roots of G _1 (ui) on the real axis or in the 
lower-half plane. In the dynamics the former correspond 
to non-decaying oscillations with frequency luq, while the 
latter are characterized by finite lifetime r = 1/r. 

Let us first consider the case when the pairing gap 



A 2 



(9) 



has a smooth maximum at the Fermi energy, i.e. the 
curvature is small < (3 « 1, and the characteristic 
energy scale 7 of the maximum is large, i.e. 7 3> A^?. 
Although in a realistic situation these two parameters can 
be related to one another, for the sake of the argument 
we consider 7 and (3 as independent parameters. In this 
case the root of G _1 (lj) = is close to the quasi-particle 
continuum w A < 2A^. By expanding Eq. J5]) in small 
parameter z = 1 — u 2 / (2A F ) 2 , we obtain 



G- 1 W«-?(^--/?ln[2 7 /A f ] 



(10) 
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FIG. 2: The analytic structure of G(u>) sketched in the com- 
plex plane of lo: as (3 changes from positive values to zero, 
the Higgs pole approaches the branch point 2Af (regime A); 
at negative values it goes into the non-physical sheet (regime 
B) through the branch cut, and first moves along the real 
axis, then obtains imaginary part at (3 < (3i\ at value (3 < (3 C 
it enters the physical sheet as a resonance (regime C). Inset 
(a): Deformation of the integration contour in Eq. (|8} with lo 
crossing the real axis, that results in the analytical structure 
shown in the main plot. 

In the dispersionless case (3 = 0, the spectral function re- 
duces to G^ 1 (z) ss — (g7r/2)i/i at z > 0, and we recover 
the square-root singularity that leads to a power-law de- 
phasing discussed in Ref. When (3 > 0, a simple pole 
of G(uj) appears on the real axis of lo at 



4Ai 1 



-J' 



In 2 [2 7 /A F 



(11) 



inside the spectral gap of quasi-particle excitations. It 
is identified as the discrete Higgs mode in regime A of 
Fig. Q] Close to the pole the spectral function 



G M ~ -j- ln [2 7 /A F ] 

n z g LO z 



(12) 



defines the spectral weight of the Higgs mode in the dy- 
namics of SA(t). As expected, it vanishes in the disper- 
sionless limit (3 = 0. The pairing dynamics at (3 > is 
characterized by non-decaying oscillations of 8A(t) with 
frequency (fTTj) . 

By inspecting Eq. (fTTj) we find that a virtual state 
is formed at smooth negative dispersion. The analyti- 
cal continuation to negative (3 < shows that the mode 
seemingly has the same energy as for (3 > 0. In fact it 
belongs to a non-physical sheet of lo (regime B in Fig. [T|) 
as it follows from careful analysis of Eq. (fT0|) . Indeed the 
decrease of (3 from positive values to zero, brings the en- 
ergy of the mode to the branch point at lo = 2Af- Upon 
further reduction of (3 to negative values, the mode enters 
the non-physical sheet of the Riemann surface parameter- 
ized by lo from the branch cut and moves away from this 



point along the real axis as a virtual state (see Fig. [2]) . In 
the pairing dynamics, the virtual state leads to a power- 
law decay of 5A(t) similar to the dispersionless case of 
Rcf. [U. 

In the opposite limit of sharp dispersion, 7 <C Ap and 
|/3 1 ^> 1, this virtual state becomes a resonance within the 
quasi-particle spectrum with finite lifetime, lo a = LOQ — iT. 
To demonstrate this effect we study a specific example of 
the energy-dependent pairing gap: 



A 2 



(3 



1 + e 2 / 7 2 



/ e = A e /A F , 



(13) 



where (3 is a curvature at the Fermi energy, and 7 is 
a characteristic energy scale. We assume \(3\^ 2 <C A F , 
so that the asymptotic value of the pairing gap A* = 
Ap{\ — /?7 2 /A 2 ? ) 1 / 2 is close to Ap. In the case of smooth 
dispersion, < (3 <C 1, we obtain Eqs. (fTT|) and (fT2|) . In 
the opposite limit of sharp energy dispersion (3 3> 1, the 
energy of the Higgs mode is close to the asymptotic value 
of the spectral gap 2A*, as it follows from the expansion 
ofEq. ©: 



4Ai 



-a; 



Jo°° de ■ 



(14) 



-A 2 -A 2 



In this expression, the upper limits of the integrals are 
substituted with infinity, by using their fast convergence. 
Explicit calculation for model (fT"3f reveals a square-root 
singularity in the energy of the Higgs mode 



4A 2 + 4 7 2 ^, 



(15) 



as a function of (3. This singularity is related to the form 
of interaction in Eq. (fl"3|) , and is not universal. 

The spectral weight of the Higgs mode in the dynamics 
is extracted from the residue of the pole in the spectral 
function 



(16) 



vanishing in the dispersionless limit 7 = (see Eq. (fT5|) ). 

The analytic structure of Eq. (fT5")) is different from 
Eq. (fTT)) : as the sign of (3 changes from positive to nega- 
tive, the former acquires an imaginary part with the real 
part simultaneously pushed in the quasi-particle contin- 
uum, since A* > Af- To verify this qualitative obser- 
vation one needs to continue G(lo) analytically in the 
lower-half of complex u-plane and locate the pole. 

The analytical continuation of G{lo) is achieved by iso- 
lating the singularities of the integrand in Eq. ([8]). We 
notice that wherever the imaginary part of frequency S lo 
changes from positive to negative value, provided its real 
part satisfies ^Rlo > 2Ap, a singularity occurs in the inte- 
grand of G _1 (w) at the energy e that solves the equation 
E e = ^Rlo/2. By deforming the contour of integration 



4 



over the energy to enclose the pole, as Fig. [2ja) illus- 
trates, we find the analytical continuation. 

Next we apply this method to calculate the energy of 
the Higgs resonance in the case of negative dispersion 
at \(3\ ^> 1. The explicit calculation of G(uj) for model 
interaction (|13[) is straightforward but tedious. Solving 
for the zeros of the analytic function, we define the energy 
of the Higgs resonance 

cul «4A 2 - l 4 7 2 ^3|. (17) 

As we have anticipated, the energy of the mode coincides 
with Eq. (|15p analytically continued to negative values 
of (3. The spectral function is obtained from Eq. (Til))) by 
substituting /3 1 / 4 with l/^/V 7 ^ 4 . 

The analytical structure of G(lo) for the energy disper- 
sion (fT3|) is shown schematically in Fig. [2l At smooth and 
negative dispersion below a critical value, |/3| < |/3 C (7)| ^ 
1, the Higgs mode belongs to the non-physical sheet 
(regime B). The critical point /3 C ( 7 ) depends weakly on 
the energy scale 7 <C Ap. In this parameter range there 
is an additional critical point |/3i(7)| < |/?c(7)| such that 
at |/3 1 < \ j3i\ the energy of the mode is real. Conversely at 
\(3\ > \/3i\ the energy develops imaginary part increasing 
until the mode enters the physical sheet at (3 = (3 C where 
it becomes a resonance (regime C). 

The change from analytic to non-analytic dependence 
on f3 in the Higgs mode energy (fTTj) and ([T5|) is accom- 
panied by a singular behavior of the quasi-particle de- 
phasing. Indeed, by expanding the quasi-particle energy 
at small e we find E e « A F + (1 - /?)e 2 /(2A F ). At 
(3 = 1 the quadratic term in e vanishes, and the disper- 
sion is controlled by higher-order terms in the expansion. 
This singularity affects the quasi-particle dephasing in 
the pairing amplitude SA qp (t) — 3?/ deu e e _l2Se *, where 
u e is a smooth function finite at the Fermi energy. Upon 
substitution of the expansion of E e to the phase factor, we 
obtain a 1/2 power-law decay SA qp (t) ~ 5fte~ ,;2AF *[(l - 
/3)t] -1 / 2 valid at long times t > 1/Ap. Function SA qp has 
a square-root singularity with a branch cut along the pos- 
itive real axis at (3 > 1. At smaller values the function is 
analytic. There is no direct correspondence between this 
non-analyticity controlled by the low-energy expansion 
and the non-analyticity of the Higgs mode (fl~5|) defined 
by the high-energy asymptotics of the interaction. The 
explicit calculation in Eqs. (fTS)) and (fl7|) illustrates that 
the character of the latter is specific to the model. 

The dynamics of the pairing amplitude 6A(t) is con- 
trolled by G(uj). The two contributions, from the Higgs 
pole w A = loq — iT located on the physical sheet, and from 
the quasi-particle branch cut, compete with each other. 
At relatively small times t > 1/Ap the quasi-particle 
states provide the main contribution leading to the al- 
gebraic decay of the pairing amplitude. At intermediate 
time scales, 1/Ap <t< 1/r, the Higgs mode dominates 
the dynamics characterized by the pairing amplitude os- 
cillating with frequency luq and decaying exponentially 



with the decrement r = 1/r (infinite for a discrete Higgs 
mode). Eventually, at even longer times <>t the quasi- 
particle contribution prevails, although having by that 
time a smaller amplitude. When the Higgs mode turns 
into a virtual state, the power of algebraic decay is in- 
creased compared to the power 1/2 discussed above. 

The analysis of the pairing dynamics within the BCS 
model (p} is restricted to time scales small compared to 
quasi-particle relaxation time T qp [l6| ■ At weak coupling 
and at small temperatures it is exponentially suppressed, 
which validates our result for the lifetime of the Higgs 
mode (dU provided 7 2 > A F /(T qpy /\j3\). 

In conclusion, the energy dispersion of pairing in- 
teraction qualitatively changes collective dynamics of 
fcrmions. The suppression of the interaction at the Fermi 
level generically leads to formation of a virtual Higgs 
state with energy on a non-physical sheet of Riemann 
surface. On the other hand, the enhancement of inter- 
action at the Fermi level makes this mode discrete with 
energy located in the spectral gap of quasi-particle ex- 
citations. Depending on the details of interaction, the 
mode may also become a resonance with finite lifetime 
within the quasi-particle continuum. 

I am indebted to L. Levitov, L. Glazman, and M. Ran- 
deria for insightful discussions of this work. 
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